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Abstract 



Vacuum solutions admitting a hypersurface-orthogonal repeated principal null direc- 
tion are an important class of 4d algebraically special spacetimes. We investigate the 
5d analogues of such solutions: vacuum spacetimes admitting a hypersurface-orthogonal 
' multiple Weyl aligned null direction (WAND). Such spacetimes fall into 4 families deter- 

• mined by the rank of the 3x3 matrix that defines the expansion and shear of the multiple 

WAND. The rank 3 and rank cases have been studied previously. We investigate the 
Q\ I 2 remaining families. We show how to define coordinates which lead to a considerable 

•n ■ simplification of the Einstein equation with cosmological constant. The rank 2 case gives 

warped product and Kaluza-Klein versions of the 4d Robinson- Trautman solutions as 
well as some new solutions. The rank 1 case gives product, or analytically continued 
Schwarzschild, spacetimes. 



X ■ 1 Introduction 

5-j 



The technique for solving the 4d Einstein equation that has received most attention is search- 
ing for solutions with algebraically special Weyl tensor [1]. This is how the Kerr solution 
was discovered. It is natural to exploit the same technique to solve the Einstein equation for 
d > A spacetime dimensions, especially as this is one of the few techniques which works with 
a cosmological constant. This paper will consider only Einstein spacetimes, i.e. solutions 
of the vacuum Einstein equation with a cosmological constant. We will look for solutions 
which are algebraically special in the classification of Coley, Milson, Pravda and Pravdova 
[2]. Such a solution admits a multiple Weyl aligned null direction (multiple WAND), the 
higher-dimensional generalization of a repeated principal null direction. 



*hsrlOOO@cam. ac.uk 
^ aahg2 @cam .ac.uk 
^cpt39@cam. ac.uk 



1 



In 4d, the Einstein equation simplifies considerably when the repeated principal null di- 
rection is hypersurface-orthogonal. There are two families of such solutions: the Robinson- 
Trautman (RT) family and the Kundt family pQ . RT spacetimes are defined by the property 
of admitting a null geodesic congruence with vanishing shear and rotation, but non-vanishing 
expansion. This family includes the Schwarzschild solution, the C-metric, and time-dependent 
spacetimes which approach Schwarzschild asymptotically. Kundt spacetimes are defined by 
the property of admitting a null geodesic congruence with vanishing expansion, rotation and 
shear. This family includes pp-waves, and the near-horizon geometry of the extreme Kerr 
black hole. 

Higher-dimensional RT spacetimes were investigated in Ref. [3]. It was found that these 
are significantly less rich than in 4d, containing certain simple time-independent generaliza- 
tions of the Schwarzschild solution but nothing that could be identified as a higher-dimensional 
analogue of the C-metric, and no time- dependent generalizations of Schwarzschild. Higher- 
dimensional Kundt spacetimes were discussed in Refs. [U [5] . 

A possible reason why the higher- dimensional RT solutions are less interesting than in 4d 
is that a multiple WAND in d > 4 dimensions need not be shear-free. Hence, unlike in 4d, the 
RT and Kundt families are not the most general solutions admitting a hypersurface orthogonal 
multiple WAND. For example, a Schwarzschild black string belongs to neither family. So it 
seems worthwhile investigating more general solutions admitting a hypersurface-orthogonal 
multiple WAND. Ref. [6] studied some general properties of such solutions. Coordinates were 
introduced analogous to those used for RT or Kundt solutions. The dependence of the metric 
on one of these coordinates was fully determined for the more special algebraic types. 

Recently, Ref. [7] investigated algebraically special solutions in 5d. In the case of a 
hypersurface-orthogonal multiple WAND, it was found that there exist two classes of space- 
time that are distinct from the RT and Kundt classes. To describe these classes, we need 
to introduce some notation. Recall that the expansion, rotation and shear of a null vector 
field i"" are defined as follows. Introduce a set of (i — 2 orthonormal spacelike vectors m", 
orthogonal to Define the {d — 2) x [d — 2) optical matrix of 

Pij = m\m]^da (1) 

The expansion, rotation and shear are the trace, antisymmetric part, and traceless symmetric 
part of piy We are interested in the case of a hypersurface orthogonal multiple WAND for 
which pij is symmetric. In this case, in 5d, Ref. [7] found that the eigenvalues of the 3x3 
matrix pij are either {a, a, a}, {a, a,0}, {a, 0,0} or {0,0,0}, where a 7^ 0. The first of these 
is the RT family and the last is Kundt family. We will investigate the two other possibilities. 
Both classes are non-empty e.g. they contain the Schwarzschild black string solution, and the 
product dS'i x S"^ respectively. 

We proceed as follows. The first step is to introduce coordinates adapted to the geometrical 
properties of these spacetimes. For RT or Kundt spacetimes, the procedure for doing this is 
well-known. However, in the cases of interest to us we have to work harder. In section |2] we 
show that these spacetimes admit foliations by null submanifolds which enable one to define 
a natural set of coordinates adapted to the structure of piy We then investigate the vacuum 
Einstein equation, allowing for a cosmological constant. 
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In section [3] we consider the case for which pij has rank 2. The general solution involves 
a function m similar to the "mass" parameter of the 4d RT solutions. For m 7^ we find 
that the solution must belong to one of the following classes (i) for any cosmological constant 
A, a warped product involving a 4d vacuum RT solution (equation (18^ ): (ii) for A = 0, the 
Kaluza-Klein oxidation of a 4d RT solution with a null Maxwell field (equation (!89|) ): (iii) 
for A = 0, a new class with metric given in equation (!94|l . for which the Einstein equation 
reduces to a pair of elliptic PDEs in 2 dimensions; (iv) for A < a similar class, with metric 
given in ( llUUp . involving the solutions of PDEs in 3 dimensions. For m = the Einstein 
equation reduces to a pair of equations identical to those governing type III (or more special) 
RT solutions of 4d Einstein-Maxwell theory (section l3.3p . 

When pij has rank 1, we determine all solutions explicitly in section HI They coincide with 
the solutions obtained in Ref. [S]. They are either product spacetimes {dS^ x or adS^, x H"^) 
or analytically continued versions of a generalized Schwarzschild solution (equation fll34p ). 
Examples of the latter are the Kaluza-Klein bubble of Ref. [9] and the anti-de Sitter soliton 
of Ref. [10]. 

Finally, we discuss a 6d case. Ref. [Tl] obtained necessary conditions on the eigenvalues of 
the optical matrix of a hypersurface-orthogonal multiple WAND in any number of dimensions. 
In 6d, there is the interesting possibility of eigenvalues {a, a, b, b} with a,b ^ 0, a ^ b. The 
coordinates of section [2] can be introduced for this case. However, the only solutions of the 
Einstein equation (with cosmological constant) in this class are conformally flat. 

2 Canonical form for metric 
2.1 Introduction of coordinates 

We will make use of the higher-dimensional generalization of the Geroch-Held-Penrose (GHP) 
formalism introduced in Ref. [12]. Our notation will follow the notation of that reference. 
Introduce a null basis {cq = i,ei = n, Cj = rrii}, i = 2, . . . ,d — 1 satisfying 

i'^ = = i ■ rrii = n ■ rrii = £ ■ n = 1 nii ■ rrij = Sij (2) 

Greek indices /i, z/, . . . refer to this basis. Latin indices a,b,c . . . are abstract indices. Introduce 
the following notation for the connection components [T3] 

i 

L^,u = "^u^fi, = VyU^, M^,u= Vi.(mi)^ (3) 

Consider an Einstein spacetime, in any number of dimensions, admitting a multiple WAND. 
Ref. [8] proved that such a spacetime must admit a geodesic multiple WAND so there is 
no loss of generality in restricting attention to geodesic multiple WANDs. We start with 
a couple of useful Lemmas. Note that we do not yet assume our multiple WAND to be 
hypersurface-orthogonal. 

Lemma 1. Let be a geodesic multiple WAND (not necessarily afiinely parameterized) 
such that (i) pij (defined in ([T])) has a block diagonal structure with blocks pjj and pa/s, i.e.. 
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Pia = Pal = 0; (ii) pij = X6ij for some function A; (iii) A is not an eigenvalue of pa/3. Then 
the distribution P = span{£, m/ (V/)} is integrable. 

Proof. Using (i), equation NPl of Ref. [12| with i = a, j = I gives 

Pa(3 Mio= pji Mjo= A Mio (4) 

a 

where we used (ii) in the final equality. Now (iii) implies that Mjo= 0, i.e., 

mj-{i- Vm„) = (5) 

Consider [i,mi]. First we have m,] = «^V6mi„-rm^Vfe4 = -miJ''Vfer-(l/2)m^Vfe£2 = 
using the geodesic equation for We also have ■[£, m/] = ma-(^-Vm/)— mQ,-(m/- V£) = 
—mi ■ {i ■ Vrua) — Pal = 0. These results imply [£, mj] G V. 

Using (i) and (ii), equation NP3 of Ref. [12] with i = a, j = J,k = K reduces to 

/3 a 

Pap M[JK]= A M[JK] (6) 

a 

and hence (iii) implies Mijk]= 0, i.e., 

■ (niK] ■ Vrria) = (7) 

Finally, consider [mj,mK]- We have i ■ [mj^rriK] = — 2m[/^ ■ (mj] ■ V^) = —2pKj = using 
(ii). We also have rrf ■ [mj^rriK] = —^miK ■ [i^j] ■ Vm^) = using the result just obtained. 
These results imply [mj.,mK\ G "C- This concludes the proof. 

Lemma 2. The basis vectors can be chosen to be parallelly transported along geodesies 
with tangent whilst preserving the properties (i), (ii), (iii) of Lemma 1. 

Proof. Let V"" be parallelly transported £"VaV^^ = with = V ■^ = Q {a. condition preserved 
by parallel transport). Then 

tdaVi = tVa iy ■ rrii) = M^o= V Moo +Vj Mjo= Mjo (8) 

where, in the final equality, we used Moo= yb'mi = —m'^i 'Vbia = (from the geodesic 

a 

equation). Using ([5]) (i.e. Mjo= 0) we have 

i-dVj = VjMjo, i-dVa = VpMpo (9) 

Hence if all Vi (all V^) initially are zero then all Vj (all V^) remain zero under parallel 
transport. Now define m- to be a set of basis vectors which coincide initially with mj and 
are parallelly transported along a geodesic with tangent C.. Using the result just derived, m'j 
remains orthogonal to and remains orthogonal to m/ (and both remain orthogonal to 
t). It follows that is a linear combination of ^ and the mj, and is a linear combination 
of ^ and the mp. Hence the new basis is related to the old one by a null rotation [12] about 
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^ (which does not change pij) and a spin [T2] which does not mix m/ with rrta and hence 
preserves the structure of pij. 

In the rest of this paper, we will impose the extra condition that the geodesic multiple 
WAND is hypersurface-orthogonal, which is equivalent to vanishing rotation: 

pm = 0- (10) 

Let pt, L = 1, . . . , N he the distinct eigenvalues of pij = p^ij). Choose the spatial basis 
vectors {m.j} = {^(/j), . . . ,m(7^)} so that spans the dh eigenspace of pij. Then p(jj) 

is diagonal. Lemma 1 implies that the distribution = {i,mj^} is integrable for all l. 
Therefore each eigenspace of pij determines a foliation by null submanifolds A/], of spacetime. 

Let us focus on the case in which N = 2, i.e., pij has two distinct eigenvalues. For d = 5 
this is the only possibility other than Robinson- Trautman or Kundt [7] (we also know that 
P2 = in this case although we shall not assume this yet). We can introduce coordinates 
adapted to these null foliations as follows. 

i is hypersurface orthogonal. Introduce a function u with du ^ so that i is orthogonal to 
surfaces of constant u. Using a boost [12] we can arrange that = —{du)a, which implies that 
i is tangent to the affinely parameterized null geodesic generators of the surfaces of constant 
u. 

Let the surfaces A/i have dimension n + 1 and introduce a chart 

(x°,...,x",F"+\...,y'^-i) (11) 

so that the surfaces A/i are surfaces of constant Y. These lie within the surfaces of constant 
u therefore u is constant over each surface A/i and hence cannot depend on X. Therefore 
u = u{Y). Hence we can define new coordinates {X , y^~^^ , . . . , y'^~'^ , u) where each y is a 
function of the Ys and the surfaces A/i are surfaces of constant y"'~^^, . . . , y'^~'^, u. Next pick 
some spacelike hypersurface S transverse to ^ and complete (y""'"^, . . . , ?/°'~^, u) to a chart 
(l/\ . . . m) on this surface. Finally, assign coordinates {u,r,y^, . . . ,y'^~'^) to the point 
affine parameter distance r from this surface along the integral curve of i (a geodesic) starting 
from {y^ , . . . y'^~'^ , u) . Since we already know that u and y^~^^ , . . . y'^~'^ are constant along 
such geodesies, this agrees with our previous definitions of these quantities. We now have a 
chart {u,r,y^, . . .y'^~'^) for which = d/dr, = —{du)a, and the surfaces Afi are surfaces 
of constant u,y"'~^^, . . . ,y'^~^- On any such surface, we have a coordinate chart defined by 
ir,y^, . . . ,2/"). 

Repeat this procedure for the foliation defined by A/2, using the same surface E. This gives 
a coordinate chart {u,r,x^, . . . ,x'^~'^) with the property that the surfaces M2 are surfaces of 
constant u,x^, . . . ,x"' (after reordering the x-coordinates). The coordinates u, r are common 
to both charts. 

Consider rj = dx^ A. . . Adx"' Ady""^^ A . . .Ady'^~'^. If 77 = then some linear combination of 
dx^, . . . , fix" equals a linear combination of dy"'~^^, . . . , dy'^^^. But then there is some 1-form 
other than {du)a = —C-a that is normal to both Mi and A/2, which is a contradiction. Hence 

Now C. = d/dr (in either chart) so i ■ dy^ = i ■ dx^ = hence i ■ r] = 0. Therefore 
i ■ {dr A rj) = f] ^ and hence dr A rj ^ 0. It follows that (r, x^, . . . , x", y"'~^^, . . . , y'^~'^) is 
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a good coordinate chart on surfaces of constant u and then (m, r,x^, . . . , x", y""*"^, . . . , y'^~^) 
is a good coordinate chart on spacetime with the property that surfaces Mi are surfaces 
of constant m, y""*"^, . . . , y*^"^ and M are surfaces of constant u,x^, . . . ,x"'. The distribu- 
tion spanned by {d/dr, d/dx^, . . . , d/dx'^} is tangent to Afi and the distribution spanned by 
{d/dr, d/dy"'~^^, . . . , d/dy'^~^} is tangent to A/2. Hence these distributions are orthogonal. 
The metric now must take the form 

ds^ = -g"du^ - 2dudr + guN {dx^ + g'^ du) (dx^ + g'^ du) 

+ gAB{dy'' + g'^du)(dy^ + g'-^du) (12) 

where indices M, N run from 1 to n and A, B run from n + 1 to d — 2. This is a special case 
of the metric discussed in Ref. p]. Henceforth we will use (as in the above Lemmas) indices 
I, J, . . . instead of Ji, Ji, . . . and a, (3, . . . instead of I2, J2, ■ ■ ■■ We will use m, tt,, . . . to denote 
either M or A. Note that the basis vectors rrii obey mj = = mf = = 0. 



2.2 Further simplifications 



We will show how the above metric can be simplified further by imposing the multiple WAND 
condition and some components of the Einstein equation. We will do this by solving some of 
the "Newman- Penrose" (NP) equations of Ref. \V2\ . 

We have chosen i"" to be tangent to affinely parameterized geodesies, and will take 
to be parallely transported along these geodesies. It follows that n"' also must be parallely 
transported, i.e., all basis vectors are parallely transported. 

In our parallely transported basis, equation NPl of Ref. [12] (with the WAND condition) 
after shifting r appropriately reduces to 



1 

Pi = 

r 



1 



P2 



r + a{u, X, y) 



(13) 



where d{u, x,y) ^ and we will allow for the possibility a = 00 corresponding to p2 = (we 
know d = 00 when d = 5). All equations below are valid in the limit a — )■ 00 unless otherwise 
noted. Here, and henceforth, an overbar denotes a quantity independent of r. 
Calculating pij gives 



Pij 



-m^mWrgr, 



(14) 



and hence 

Taking = m in i ■ Vmf 

mf{u,r,x,y) = 
Taking fj, = r gives 



^drgr. 



and using ([T5|) gives an equation which can be solved to give 



-mj'iu, X, y), m'^iu, r, x, y) = , 
r 1 + r/a 



drg'"" = 2{drml)m'P 



-m'^{u,x,y) 



(16) 
(17) 
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We now have 



^l^IN 1 ^M-N AB 1 ^A^B 

g = —mj rrij , g - 



^2-"v -1 . ^ (1 +r/a)2 " " ^ ^ 

and hence 

guN = r'^gMN{u,x,y), gAB = {l + r/aYgABiu,x,y) (19) 
Since di = we have 

= L,i = r,- (20) 

Equation NP3 of Ref. [12] is 

^[jP\i\k] = 0. (21) 

Setting j = J,i = I,k = K gives 

1 



r 



-m^ + Tj = if n > 1 (22) 



(recall that n is the number of x-coordinates, the restriction n > 1 arises because ( |2T]) is 
antisymmetric on jk). Setting j = J,i = I,k = j gives 



Mij = (1 + r/a) ( + ) 6ij (23) 



Setting j = I,i = a,k = (3 gives 



— (^mj + mfdMoj + Ti 



r + a 



(24) 



Ma/3 = - 

j = (3,I = a,K = 'y gives 

7 ~~T f"^a + ^a^^Ao) + = ^ Hd — 2 — n>l, (Ij^OO (25) 

(r + a) ^ ^ 

The remaining components of (12 ip are trivial. 

Since the basis is parallelly propagated, = NiQ = 0. Equation NP2 of Ref. [12] (with 
vanishing curvature terms) reduces to 

n = -fj{u,x,y), r» = / fo,{u,x,y) (26) 
r 1 + r/a 

The orthogonality relations of the basis vectors give 

Ur = 1, nm = -mlmim (27) 
Using this we can calulate Xj = Ln directly to obtain 



Ti = -drm\ + -mimrrf^mjndrg'^'^ (2^ 



Since we know the r-dependence of the LHS we now can obtain (using (fT8l) ) 

Tjfj = - fi + ^rrfi (29) 
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Now using fl22|) gives 

?Tij = if n > 1 

and (|25|) gives 

mj^ = — afa — m;^(9Aa if d — 2 — n > 1, a ^ oo 
A calculation from the definition of the LHS gives 

1 „\ . 1 



Mij = , ,^ ( -Ta + -m'a ) Sij - f .. rniMmjNm^dAg^ 
l + r/aV r / 2(1 + r/a) 

Comparing with fl23|) gives 
and hence 

gMN{u,x,y) = P{u,x,yy^hMN{u,x) 

where 

OaP « / 1 



-r V a y 

Combining with ( l32l) gives (absorbing a function of ti, X into h]\,j]\[) 

P = a lid — 2 — n>l, a^oo 

Now we calculate from the definition of the LHS 

Ma/3 = f-T/ + -mj) - ■^m^AmpB'mfdMg'^'^ 

r + a \ r J 2r 



Comparing with gives 



2 

^M^AB = --9AB (gMNrnfrfi} 
a ^ 



and hence 

gAsiu, X, y) = Q{u, x, yY'^hABi.u, y) 

where 

dhlQ 1 _ - N - r 

Q a 

Combining with f l3T|) gives (absorbing a function of u, y into /i^s) 

(5(w, x, y) = 1 if a = oo or n > 1 
Now we return to f|T7j) . This gives 

rM - r - M i -rM 

= —m^mf +g 



so, from fl3T|) . 

gvM ^ -vM i{n>l (44) 

and, from (14T|) . 



grM ^ a_ M^OnQ ^ _,M if ^ ^ (45) 



From (fTTl) we also obtain 



= n ^"/-^2 f^" + -<) + if a ^ oo (46) 

(1 + r/ay \ a J 

g"-"^ = -2rf^mi + if a = oo (47) 

Note that the coefficient in the r-dependent piece involves (in both cases) 



f„ + ^< < = -g^'^dBP/P (48) 
a J 

In summary, we have shown that the metric components in f|T2|) must take the following form: 

9mn{u, r, X, y) = PP{u, x, yy^huNiu, x) (49) 
where P = ai{d — 2 — n>l and a ^ oo 

QAsiu, r, X, y) = {1 + r/a{u, x, y)fQ{u, x, yy^hAsiu, y) (50) 
where Q = 1 if a = oo or n > 1, 

/^^(n,r,x,i/) = lf^^(«,a;,y) + r*^(n,x,i/), (51) 

where /^*^ = if n > 1 and 

JrM ^ ^p2j^MNd^ if a ^ OO (52) 

Q 

^jrM ig unconstrained if n = 1, a = oo.) 

9'^ = -7T-^j^Q'h^''^ + f^in,x,y) if Moo (53) 

(7^^ = 2rg2/i^^^ + ^-^(^x, X, y) if a = oo (54) 

where h^^^ , h^^ are the inverses of Hmn, hAB- 

In deriving these results, we used only equations NPl, NP2 and NP3 with vanishing (Ricci 
and Weyl) curvature components (as is the case for an Einstein spacetime with a multiple 
WAND). The metric deduced above is the most general metric admitting a hypersurface- 
orthogonal geodesic multiple WAND, for which p^j has precisely two distinct eigenvalues, and 

satisfies Rrr = RrM = RrA = 0. 
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2.3 The 5d case 



In 5d, Ref. [7] proved that a = oo so there are just two cases to consider, n = 2 and n = 1. 

For n = 2 we can perform a coordinate transformation x — )■ x'{x,u) to arrange that 
huNiu^x) = 6mn (absorbing a conformal factor into P). Hab has a single component h{u,y) 
and y — )• y'{y, u) can be used to set h = 1. The metric is 



^2 



ds^ = -Udu^ - 2dudr + — (dx + Adu) + 

pi 



dy+ lc + 2r^] du 



(55) 



where we have dropped the bar on P. Here bold font denotes objects transforming as vectors 
in 2d Euchdean space with coordinates x. The functional dependence is U = U{u,r,:K,y), 
P = P{u, X, y), A = A{u, x,y),C = C{u, x, y). 

For n = 1, hMN{u,x) has a single component which we can set equal to 1 using x — )■ 
x'{x, u). Although we have Q = 1 m. this case, a transformation y — y'{y, u) can be used to 
achieve hAsiu^y) = Q{u,y)~'^6AB for some function Q. The metric is 



,2 



ds = —Udu — 2dudr H - 

p2 



/A \ l'^lr / V P\ 

2 o^.,^. , dx+ — + B]du ' . / . o.^2 V \ 



1 



dy+ [C + 2rQ'^j 



du 



2 



(56) 



Again we have dropped the bar on P. Here bold font denotes objects transforming as vectors 
in 2d Euclidean space with coordinates y and V is the gradient operator in this space. The 
functional dependence is U = U{u,r,x,y), P = P{u,x,y), A = A{u,x,y), B = B{u,x,y) 
C = C{u,x,y), Q = Q{u,y). 

As mentioned above, these metrics satisfy Rrr = Rrx = Rry = 0. We will determine the 
conditions arising from the remaining components of the Einstein equation with cosmological 
constant: 

Rab = ^9ab (57) 

The calculations were performed using Mathematica. 

3 5d metric with rank 2 optical matrix 
3.1 The Einstein equation 

For the metric (l55|l . the components of the Einstein equation in the x^ ,y directions give 

Rf^^y = ^ dyA^^ = (58) 

so A = A(m, x) 

R^2 = ^ diA^ + d2A^ = (59) 
- i?22 = ^ diA^ - d2A^ = (60) 
These equations are equivalent to (lowering indices with 6mn) 

di^M^N) = ■ ASmn (61) 
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i.e. A is a conformal Killing vector field of 5mn- This implies that A can be eliminated by 
a change of coordinates x — )■ x'{x^u) and absorbing a function of u,x into P. So henceforth 
A = 0. 

We find that 

Hence equating powers of r in Ryy = A gives 

dyC = C = C(m,x) (63) 

Pyy + = (64) 

where we introduce the notation Py = dyP, Pyy = dyP etc. 

The equation Ru + R22 = 2Ar^/P^ involves U and drU. It can be integrated with respect 
to r to give (using (!64l) ) 

f/ = -— + A log P - 2r [(log P)„ - C(log P)y\ + r^'"" " 

where m = m('u, x, ?/) is the "constant" of integration and A = P'^5^'^^ dudN is the Laplacian 
of Hmn = P^'^^MN- Note that the r-dependence of the metric is now fully determined. 

Using these results, we find that the ur component of the Einstein equation is satisfied. 
Next we consider the uM components. These have terms that are (9(l/r^) and 0{r). The 
former give 

c^A/m = ^ m = m{u, y) (66) 

The 0{r) terms vanish as a consequence of flM|) . The uy component of the Einstein equation 
has terms that are (9(l/r^), 0{1) and 0{r). Using f l64|) . only the (9(l/r^) part is non-trivial, 
giving 

AC = -2my + 8m(log P)y (67) 

If we divide through by P^ then the LHS is independent of y. Taking a y-derivative and using 
(!64|) we obtain 

12m(log Pfy - 6m3^(log P)y + rUyy + Am = Q (68) 

The final component of the Einstein equation is the uu component. Using the above solution 
for U and equations ( l64l) . ( !67l) . ( l68l) . we find that this reduces to 

AAlog P + 12m [(log P)u - C(log P)y] - A{mu - Cruy) - P'^iVCf = (69) 

where {VCf = d^'^duCdNC. 

In summary, the Einstein equation reduces to equations ( l64|) . ( 1671) and ( l69l) with C = 
C{u, x) and m = m{u, y). 

Before analyzing these further, we give the expression for the boost- weight zero Weyl 
components [12]. Choosing the spatial basis vectors m2 = {P/r)d/dx^, = {P/r)d/dx'^ 
and m4 = d/dy gives 

—m/r^ 

$ij =1 -m/r^ I (70) 

P9iC/2r2 Pd2C/2r^ 
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From this we can read off $ = —2m/r^ and ^fj^fj = P^(VC)^/(8r^). These are invariants of 
the solution, i.e., independent of how i is completed to a basis. vanishes (which implies 
that the solution is type III or more special) if, and only if, these two scalar invariants vanish. 
If 7^ then there is a curvature singularity at r = 0. 



Fyy +-F = 0. (71) 



3.2 mj^O 

In this case, (|68l) is a quadratic equation for (logP)^ whose coefficients depend only on u,y. 
Hence (logP)y is independent of x. It follows that P has the form P = F{u,y)G{u,x.). 
Equation flM|) gives 

A 

4 

This is trivial to solve. Using the gauge freedom y ^ y + f{u), C ^ C — f'{u) we can bring 
any solution of this equation to the form a{u)F{y). The function a can be absorbed into G. 
Hence, without loss of generality (and dropping the tilde) 

P = F{y)G{u,^). (72) 

where F{y) obeys (ITTll . Plugging into equation (J6711 now gives 

F^(|^)^ = -1g^V^C (73) 

where = S^^^dudN (so A = P^V^). The LHS depends only on -u,?/, the RHS depends 
only on M, a; hence we must have 

v^C = -]{u) (74) 

for some function j{u). Integrating now gives 

m(w, y) = ]^F{yY {j{u)H{y) + fc(n)) (76) 

where 

Hiy) = J F{y)-^dy (77) 
Now consider equation fl69|) . Dividing by P^ gives 

G'V' (G^V^ log g) + {jH + k) (q^ + 2G^^ - 2 {j^H + K) + 2^ - ^(VC)^ = (78) 

Now differentiate with respect to y, use flTTj) . and multiply by P^ to obtain 

= J (e^^ - 2G^) - 2G{jH + k) [fI + ^F') - 2j. + 2^(VC)^ (79) 
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Differentiate again with respect to y, using fl7T]) to obtain 

= -2^ [f^^ + ^F') {VCy (80) 

There are two cases to consider. 

Case 1. VC = (equivalently $^ = 0) hence C = C{u). This is the only case if A > 0. 
Equation fITT]) imphes that 

+ JP'' = l (81) 
for some constant A. From ( FM|) we have i{u) = 0. Equation f l79l) reduces to 

= -'^CkF-^X (82) 

hence C = (A; 7^ since we're assuming m 7^ here). Equation (ffSll reduces to 

G^V^ (g^ log g) + 6A;(log G)„ - 2A;„ = (83) 

This is the 4d Robinson- Trautman equation with "mass function" k/2. If we now perform a 
coordinate transformation r = F{iiYR then the metric becomes 

ds^ = dy^ + Fiyf \-Viu, R, x.)du'^ - 2dudR + — ^— rrdx^ ) (84) 

V g(m,x)^ ; 

where 

= + GV^logG - 2i?(logG)„ - ^ (85) 

The 5d metric is a warped product of a 4d Robinson- Trautman solution with cosmological 
constant A with the y-direction. The warp factor F{y) is obtained from flSTl) . If A > then 
we must have A > and, shifting y hj a constant, the solution is 

Fiy) = ^sm{VXy/2) (86) 

It is easy to understand why this is a solution: we can write the 5d de Sitter metric as 

ds^ = dy^ + Fijjfdsl (87) 

where ds\ is the 4d de Sitter metric with cosmological constant A. The metric ds\ can be 
replaced with any other 4d Einstein spacetime with the same cosmological constant. Our 
solution corresponds to choosing this 4d spacetime to be a Robinson- Trautman spacetime 
with cosmological constant. 

For A = we have A = or A > 0. In the former case, F is constant so we can take 
F = 1 and the solution is simply the product of a 4d Ricci-flat RT solution with a flat 5th 
direction. In the latter case we have F{y) = yjx/3y (shifting ?/ by a constant). In this case, 
the 5d solution corresponds to writing Minkowski space in the form ( IHTl) with dsl the 4d de 
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Sitter metric with cosmological constant A, then replacing dsl by a 4d RT spacetime with 
cosmological constant A. 

If A < then we have have A>0,A = OorA<0 and the solutions are obtained by 
writing 5d anti-de Sitter space as a warped product ( 187|) where dsl is 4d de Sitter, Minkowski 
or anti-de Sitter space, then replacing dsl "with a 4d RT solution with the same cosmological 
constant. 

Case 2. In this case, 

+ (A/4)F2 = (88) 

This is possible only if A < 0. We consider A = and A < separately. 

Case 2.1. A = 0. We can set F = 1, (absorbing a constant into G) and then H = y. 
There are two subcases to consider depending on whether or not j = 0. 

Case 2.1.1. j = 0. Then m = k{u)/2 and the metric is 

ds"^ = -U(u, r, x.)du^ - 2dudr + — — rfx^ + (dy + C(u, x.)duf (89) 

G{u, x)^ 

where 

f/ = -M!^ + G'V'logG-2r(logG), (90) 

Note that d/dy is a Killing vector field in (189|) and this metric takes the form appropriate to 
a Kaluza-Klein reduction to 4d. The KK scalar field is constant and the KK Maxwell field is 
obtained from the 1-form potential C{u, x)du. This is a null field, as it must be for consistency 
with the scalar field equation of motion. Equation (fMll gives the 4d Maxwell equation: 

V^C = 0, (91) 

and the 4d Einstein equation is ( 1781) 

[G^VHogC) + 6A;(logG)„ - 2A;„ = ^^(VC)^ (92) 

This is the 4d RT eq with a null Maxwell field. So this 5d solution is simply the KK lift of a 
4d RT solution of Einstein-Maxwell theory, with a null Maxwell field. From fITOl) we see that 
this is an example of a solution with a hyper surf ace-orthogonal multiple WAND and 7^ 0. 

Case 2.1.2. j ^ 0. Note that this implies VC 7^ and hence $^ 7^ 0. We can use the 
gauge symmetry u — > u{u'), r — > r/du'U to rescale m by an arbitrary function of u. This can 
be used to set j = constant. (We could set j = I but it is convenient to allow more general 
values e.g. so that the limit j — )■ can be understood.) Equation (179ll reduces to G„ = and 
hence 

a = G(x) (93) 

In the analysis above, we fixed the gauge symmetry y — )■ y + f{u), C — )■ C — f'{u) by imposing 
fl72]) . However, in the present case of constant F, this gauge symmetry remains and we can 
use it to set = (since j ^ 0) so m = y /2. The metric is 

ds^ = -U{r, X, y)du^ - 2dudr + — — rfx^ + {dy + C{u, x)duf (94) 

G(x)^ 
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where 

U = -— + G^y^\ogG (95) 
r 

The Einstein equation reduces to ( !74|) (with j = constant) and (178!) . The latter simphfies to 

G'V' (g'V' log G) = G\VC)^ - 2jG (96) 

As with all the m ^ solutions, this solution cannot be type III or more special (since 
$ij 7^ 0). It cannot be type D either: it can be shown that it is impossible to achieve = 
by a null rotation about i"-. Since G is independent of u, the above equations might force G 
to be independent of u. In this case, d/du is a Killing vector field and the metric also admits 
a homothety: rescaling r,y,u by the same constant rescales the metric by a constant. 

Case 2.2. A < 0. Write A = -A/L"^, L > 0. Equation ([HHD gives F{y) = e^/^ where we 
have absorbed a constant coefficient into G and used the freedom y — )■ —y, (with G — t- —G, 
j -j) to arrange F' > 0. We can take H{y) = -(L/2)e"^^/'^. Equation ([79]) is 



G^ 

Using this, equation (178|) becomes 



2_\ ^Ir 

GVu L 



G^VGy-jG (97) 



G^V^ (G^V^ log g) + 6fc(log G)„ - 2K = -jkG (98) 

If j = then (1971) implies G = G{u), which is case 1. So we must have j ^ here. Under the 
gauge freedom y ^ y + f{u), G G — f'{u), G — )■ Ge^^^^ we have 

J ^ e^^/^j, k ^ e^^/^A; (99) 

We also have the gauge freedom u — )■ u{u'), r — )■ r/du'U which rescales m by an arbitrary 
function of u and hence rescales j, k by the same arbitrary function. These transformations 
can be used e.g. to set j = I and k = constant. 

A coordinate transformation r = F{y)'^R brings the metric to the form 

ds'^ = (dy + G(u, y.)duf + e^^/^ \-Viu, R, x, y)du^ - 2dudR + — rfx^ ) (100) 

V G(m,x)2 / 



where 



V = -^( k{u) - Mlie-^yl^\ + G^V^ log G-2R [(log G), + ^ 



R \ ^ ' 2 



L 



(101) 



The Einstein equation has been reduced to the three equations ( 17il) . f l97|) and ( l98i) . 

The above metric appears to be asymptotically locally anti-de Sitter as ?/ — )■ oo. The 
metric on the conformal boundary is similar to a 4d RT solution of Einstein-Maxwell theory 
(without cosmological constant) where the Maxwell potential is Gdu. However equation (|74l) 
shows that j plays the role of a source for G . Maybe one can interpret the boundary metric 
in terms of a RT solution of Einstein-Maxwell theory coupled to null charged matter. 
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3.3 m = 



Recall that m = is a necessary, although not sufficient, condition for the spacetime to be of 
type III or more special. When m = 0, the Einstein equation reduces to flM|) along with 

AC = (102) 

AAlog P = P\VCf (103) 

These two equations are identical to the equations governing d = 4 type III (or more special) 
RT solutions of Einstein-Maxwell theory [1] (with the Maxwell field obtained from the poten- 
tial Cdu). The difference here is that P has y- dependence, determined by f l64l) . The latter 
equation can, of course, be solved explicitly e.g. P = q;(m,x) -|- I3{u,:s.)y when A = 0. Recall 
that C is independent of y. 

From ( !70|) . the (5d) solution is type III or more special if, and only if, VC = 0. In this 
case, C = C{u). Using the gauge freedom y ^ y + f{u), C —> C — f'iu) so we can set C = 0. 
In this case the above pair of equations reduces to AAlogP = 0, which is the same as the 
equation governing d = A type III (or more special) vacuum RT solutions [Tj. Again, there is 
a difference in that P has y-dependence, restricted by (1641) . 

With C = 0, a calculation of "^[jk reveals that the solution is type N (or O) if, and only 
if, we have 

V(AlogP)=0 ^ A\ogP = K{u,y) (104) 
for some function K and 

= (105) 




where we have defined a complex coordinate C, = {x^ + ix^)/\/2 and d = d/dC,. Equation 
f ll04p is the same condition as arises for a 4d vacuum RT solution to be type N (or O) [1]. It 
implies that 2-surfaces of constant m, r, y have constant curvature. 

Let us return to the most general m = case and attempt to solve the equations by an 
Ansatz motivated by our results for m 7^ 0: 

P{u,^,y) = F{u,y)GM (106) 

so (IMI) becomes (fTTj) . Equation (11031) implies either dyF = or VC = 0. In the first case we 
must have A = and the solution is the same as case 2.1.1 above but with m = (fc = 0), 
i.e., the KK lift of a 4d RT solution of Einstein-Maxwell theory with a null Maxwell field. In 
the second case the solution is type III (or more special) and, as explained above, we could 
use the freedom y ^ y + f{u) to set C = 0. However, we will keep C{u) 7^ and instead use 
this freedom as we did below (I7T]1 to arrange that F = F{y) where F satisfies f lHTj) for some 
constant A, as in case 1 above. A change of coordinates r = F^R then brings the metric to 
the form 

ds'^ = (dy + C(u)duf + Fiyf \-Viu, R, x, y)du^ - 2dudR + — 5— ^rc/xM (107) 

V G(m,x)2 ; 
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where 

V = G^VHogG - 2R[{\ogG)u + C(logF),] - ^ (108) 
and the single equation of motion is 

(C^V^logG) = (109) 

When C{u) = this solution corresponds to setting m = (fc = 0) in the warped product 
solutions of case 1 above. 

4 5d metric with rank 1 optical matrix 
4.1 The Einstein equation 

Now consider the metric fl56|) . As mentioned above, the rr, rx and ry components of the 
Einstein equation already are satisfied. The xy^ components are 

3r 1 / 1 4A 1 \ 

= = -i^^doB + - [-p-,dnA + —dnP + g^d^Cn^ (110) 

The 0{r) term gives OdB = hence B = B{u,x). But this implies that we can use a change 
of coordinate x = x{x',u) to set 

B = 0. (Ill) 
The 0{l/r) terms give {Cb = SbaC^ = C^) 

O^Cb = P'Q'Ob (-^) (112) 

The AB components of the Einstein equation can be treated similarly. They are equivalent 
to the equations 

diACB) = -^D^QSab (113) 



d^A {Q^dB)P) = \qvp ■ vq + ^ (g' log g - A 



Sab (114) 



3 

where V is the gradient with respect to y^ and ■ denotes the scalar product defined by 6ab- 
We have also introduced a modified derivative: 

D^ = d^~C-V (115) 

These equations show that C and Q^VP are conformal Killing vector fields of Sab- However, 
we cannot remove C by a coordinate transformation of the form y = y{y',u) because C 
depends on x. 

The trace of flll4p gives 

V^P 2 /_2, „ A 



V^logg-— (116) 



P 3 \^ ■ ^ Q2 
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Next we consider the xx component of the Einstein equation. This involves U and Ur with 
the r-dependence of all other terms explicit. Integrating gives 

U = ^,+m- 2rD. log P + {^^^{VPf - gv^P - ^) (117) 

where the arbitrary function m = m(-u,x,y) plays the role of the "constant" of integration 

Using the above results, we find that the ur component of the Einstein equation is satisfied, 
and the Ricci scalar obeys P = 5A. If we let Eat = Rab ~ ^9 at then the contracted Bianchi 
identity now implies 

dr {rE^,) = dr {tE^a) = (118) 

and hence only the 0{l/r) terms in the ux and uA components of the Einstein equation can 
contain new information. The uA component reduces to an equation for which can be 
integrated to give 

m = M(m, x)p2 - {d,^Pf + 2PdlP (119) 
for some function M(m,x). Using this in the ux component gives 

2PSilP + 9, (MP2) = -2PDu (120) 

The remaining component is the uu component. The contracted Bianchi identity reveals that 
only the 0{l/r) terms in this can be non-trivial. These terms give: 

duM = ... (121) 

where the ellipsis denotes a lengthy expression involving P, Q, A and their derivatives. 

It is convenient to introduce a complex coordinate C = (?/^ + ^l/^)/\/2 so that dy'^ = 2d(d(. 
The (( component of equation flllSp gives 

BC^ = (122) 

where d = d(^, and hence 

C'^ = W{u,xX) (123) 
for some holomorphic function W. The (( component of the same equation gives 

duQ^^ = d {Q'^w) + B {Q-^w) (124) 

Differentiating with respect to x gives 

d (Q'^W^ + B (Q-^W^ = (125) 



We also have f lll2p which gives 

= P"Q"9 (^) (126) 



Note that the term in U cancels with another term contributing to g^u- This imphes that the non- 
trivial metric components are quadratic or linear functions of r, in agreement with Ref. 6 . This reference 
assumed — 0, which is true for the case considered in this section ^. 



18 



together with the complex conjugate of this equation. 
The CC component of equation flll4p gives 



d (Q^dP 



(127) 
and hence 

Q^dP = 7 (128) 
for some function J-'(u, x, C) depending holomorphically on C,. 

4.2 Invariants 

From Ref. [7], we know that there exists a basis for which $jj = diag(— $, $, $). We find 
that 

^ = ^^-\Q^dd\ogQ (129) 

If $ = then = so the spacetime must be type III, N or O. But a spacetime of type N 
or (in 5d) type III cannot have a rank one optical matrix [13]. Hence $ = if, and only if, the 
solution is conformally fiat, i.e., locally isometric to (anti-)de Sitter or Minkowski spacetime. 

Now assume $ 7^ and consider the boost weight —1 Weyl components and ^'ij]^- Use 
the basis 

P d d d 

ms = — mg = g— , m4 = Q-^ (130) 
r ox oQ oQ 

It can be shown that the equations derived above (in particular fll26p ) imply that = and 

% = Kp^Qd[p-^<^) (131) 

with \I'4 obtained by complex conjugation. Consider a null rotation about i"' with parameters 
Zi. This leaves invariant (and hence zero) but \E'3 — t- ^3 — 2$Z3. Hence we can eliminate 
\E'3 by an appropriate choice of z^. Its conjugate \I'4 also will vanish. So a null rotation with 
arbitrary Z2 can be used to set \E'' = 0. A calculation using the above equations then reveals 
that, in the new basis, "^ijk = 0. Hence this null rotation has eliminated all of the boost 
weight —1 Weyl components. 

Next we compute Q'^j in this new basis. The result is independent of the arbitrary param- 
eter Z2 that specifies the basis. Using the above equations (including fll2ip ) we find that the 
only non-vanishing components are 

^^'23 = 2$ [rP-^Q [PdP^ - P^dP) + PQ-^W^^ (132) 

and components related to this by complex conjugation (1^24) or symmetry. The solution is 
type D if, and only if, $ 7^ and VL\- = vp^^ = 0. Hence 

$ ^ 0, PdP^ - P^dP = 0, W^ = ^ type D (133) 

If the type D condition is satisfied then we can employ an argument of Ref [7] . We have noted 
that the type D condition is satisfied in a basis which, at any point, is specified by an arbitrary 
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parameter Z2- We can choose Z2 to be an arbitrary function of spacetime position. Since the 
type D condition is satisfied in this basis, n°- must be a muhiple WAND. But the vector 
field depends on the arbitrary function Z2. n°- will be geodesic only for special choices of 
function 2:2, in general it will be non-geodesic. Hence our spacetime admits a non-geodesic 
multiple WAND. All 5d Einstein spacetimes admitting a non-geodesic multiple WAND were 
determined in Ref. [8]. The possibilities are: dSj, x S*^, AdSj, x and the metric 

ds" = r'dSl + ^ + f{r)dz\ f{r) = k - - (134) 

where m 7^ 0, A; G {1, 0, —1}, ds^ is the metric of a 3d Lorentzian space of constant curvature 
(i.e. 3d Minkowski or (anti-)de Sitter) with Ricci scalar 6k and the coordinate r takes values 
so that / > 0. This is an analjd^ically continued version of the Schwarzschild metric (allowing 
for a cosmological constant). 



4.3 Solving the Einstein equation 

We will show that the only solutions in this class are either conformally flat, or the type D 
solutions just discussed. 
Case 1. ddP = so 

P = g + g (135) 
where Q{u,x,() is holomorphic. Then we have 

gV^iogg = A (136) 

This is Liouville's equation. It says that surfaces of constant u, r, x have constant curvature 
A. 

Case 1.1. A = 0. Using the above equation in (11291) gives $ = so the solution is flat. 
Case 1.2. A 7^ 0. The general solution of Liouville's equation is 

^5 (l + AZZ/2)^ , , 

for some function Z{u, (). A transformation ( — )■ Z{(, u) can be used to bring Q to the form 

Q = l + ACC/2 (138) 

Using this, equation f ll27p gives 

(1 + ACC/2) d'^G + KCdG = (139) 

which implies dQ = (e.g. act first with d) hence dP = so P is independent of ( (and (). 
Now flT26|) gives 

Q-^W^^ = d{A/P^) (140) 

and plugging this and its complex conjugate into fll25p gives dd{A/P^) = hence A/P^ = 
T-L + T-L for some T-L{u, x, (). f ll40p then gives Wx = Q^dl-L. Taking the complex conjugate and 
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acting with d gives d{Q^d'H) = 0. This imphes &H = using the same argument we used to 
show dQ = 0. Complex conjugation gives &H = and hence = 0. 

These results imply that the solution satisfies the type D condition derived above. Hence 
it is one of the solutions of Ref. [8]. From the fact that it has constant $ it follows that it 
must be dSs x ^Mf A > or AdS^ x ii A < 0. 
Case 2. ddP ^ 0. Let (, () be real and satisfy 

ddx = (141) 

We then have dP = Tddx = d{Tdx) and hence 

P = ^dx + H (142) 

for some Ti^UjxX)- Equation (11161) is 

ddP = -i (9(91ogg-2 + Ag^2) P (143) 

Using the above expression for P, multiplying by jT and using the definition of x gives 

- i-L dJ- 

d log ddx -qdx = qy- ^ (144) 

where 

g(n, C, = (Q^dd log + A) (145) 
Differentiate with respect to x: x and Q are independent of x so we obtain 

Now act with d to obtain 

{dq)(^ =0 (147) 

Case 2.1. 9g = hence (as q is real) dq = Q. Act with d on fll44p and recall the definition 
of X to obtain 

dd\ogQ-^ - qQ-"^ = Q (148) 
Plugging in the definition of q gives 

8Q^aaiogQ-A = ^ $ = (149) 

Hence this case is conformally fiat. 

Case 2.2. {'H/T)x = 0. Hence "H = fiT for some function h{u,()- From fll46p we have 
{d log F)x = 0, which implies = ^o{u, x)Fi{u, Q. We now have P = (j)x + ^) • Bnt note 
that the definition of x determines x only up to addition of a real harmonic function hence 
we are free to make the redefinition X ^ X ~ 9 ~ 9 for an arbitrary function g{u, Q. We can 
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use this freedom to eliminate h. So Ti = and P = Tdx. Reality of P implies that arg J-'q is 
independent of x. By absorbing a phase into T\ we can assume that J-q is real. 

Now consider the effect of a holomorphic change of coordinates C = Z{u^ ('). In the primed 
coordinates we have J-"' = T jd' Z. Hence by an appropriate choice of Z we can arrange that 
T'x = \. Henceforth we drop the primes and subscripts, and so 

J= = 7 = 7{u,x), P = J=dx (150) 

Reality of P implies that 

dx = Bx (151) 

Hence Q, P are independent of Im(^. Using this, fll25p can be rewritten as 

{W^ + W^)d\ogQ-'^ + dW^ + dW^ = Q (152) 

Act with dd and use holomorphicity of W to obtain 

{W^ + W^.)5' logQ-' + [dW, + dW^) logQ-2 = (153) 

We can view this pair of equations as simultaneous linear equations for Wx + and dWx + 
dWx. 

Case 2.2.1. This linear system has a non-trivial solution. This requires vanishing determi- 
nant: 

(d log (d^ log Q-^) - log = (154) 
This can be integrated to give 

log Q-^ = t{u)Q-^ (155) 

for some function t, which implies 

d'^\ogd'^X = td'^X (156) 

Integrating gives 

d^X = e*^^' (157) 

for some function s{u). Using the above results, equation fll44p reduces to t = q. However, 
from the definition of q we have q = —{t + A)/3. Hence q = —A/4 which implies $ = so 
this case is conformally fiat. 

Case 2.2.2. The linear system f ll52p . fll53p admits only the trivial solution, so Wx + Wx = 0, 
which implies 

Wx = ia{u,x) (158) 
for some real function a. Equation fll26p gives 

Hpl) = P^ (159) 

which implies 

dd (4l) = (P-'Q-') (160) 



VpV 
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But note that ad {P ^) is real so reality of the LHS implies 



ad[P^^Q-^)=0 (161) 
Case 2.2.2.1. d {P^'^Q^'^) = 0. Writing P,Q in terms of x gives d'^{dx)^^ = hence 

where 7(m) and S{u) are real. One can now determine P,Q and plug into fll44p . The final 
result gives $ = 0, i.e., a conformally flat solution. 

Case 2.2.2.2. a = so Wx = 0. Together with f llSOp . this implies that the type D condition 
is satisfied if $ 7^ 0. So in this case, the solution is either type D or conformally flat. As 



explained above, if it is type D then it is one of the solutions found in Ref. [S], i.e., dS^ x S^, 
adS^ X H'^ or the solution fll34p . 



5 A 6d case 

Ref. [llj obtained necessary conditions on the eigenvalues of the optical matrix of a hypersurface- 
orthogonal multiple WAND in any number of dimensions. We will investigate the possibility, 
in 6d, of eigenvalues {a, a, b, b} with a,b 0, a b. 

In the notation of section |2l take d = 6, n = 2 and finite a (i.e. pij has two distinct pairs of 
coincident non-zero eigenvalues). Using coordinate transformations x — )■ x'{x, u), y ^ y'iv, u) 
we can arrange that 

hMN = P{u,x)~^Smn hAB = Q{u,y)~'^SAB (163) 

Note that P and Q are not the same as P = a and Q = 1 (of section |2]). The metric can be 
written 

ds"^ = —g^'^du'^ — 2dudr + — — (dx + A(u, x, y)du) (164) 

a[u, X, yYP[u^ xY 

dy + {C{u, X, y) — Q'^{1 + r/d)~^Vyd^ du 

where we are using 2d vector notation x = x^ , A = A'^'^ , VyO = etc. Using Mathematica 
we now calculate the Ricci tensor via the Ricci identity. We find that the rr, rx and ry 
components vanish. Vanishing of the xy components gives 

The traceless part of the xx components gives 

Bi^mDn) cc 6mn, d(^M (-P^f^7V)«) oc 6mn (166) 

where Dm = D^^ . Hence D^'^ and P'^dMd are conformal Killing vector fields of 5a/ at. Similarly 
the traceless part of the yy components reveals that and Q'^d^d are conformal Killing 



+ 1 + ^ Q{u,y) 
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vector fields of 5ab- We can now eliminate by a change of coordinates = x^\x',u) 
and we can eliminate E"^ by y"^ = y"^{y',u). So D'^^ = E"^ = henceforth. 

Next we consider the trace of the xx and yy components of the Ricci tensor (assuming 
A = now). These are both linear in g^^ and drQ^^ hence we can solve to determine g^^ and 
drg^^. Consistency of the resulting expressions gives the following equations: 

R. + Ry = 0, (167) 

where Rx is the Ricci scalar of P^'^Smn and Ry the Ricci scalar of Q'^Sab- This implies 
that Rx and Ry are both functions only of u hence we have 2-spaces of constant curvature (of 
opposite sign). We also have 

2^ = Q'Vja + aRy (169) 

Calculating the Weyl tensor using Mathematica and making use of the above equations, we 
find that the boost weight zero Weyl components all vanish so the solution is type III or more 
special. However, a type III or N solution with hypersurface orthogonal multiple WAND 
cannot have a rank 4 optical matrix [13]. It follows that the solution must be conformally 
flat. 
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